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Abstract. We describe the hadronization of quark matter assuming that quarks 
C^l ' creating hadrons coalesce from a continuous mass distribution. The pion and 

antiproton spectrum as well as the momentum dependence of the antiproton to pion 
\q \ ratio are calculated. This model reproduces fairly well the experimental data at RHIC 

■ energies. 
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1. Introduction 

X 

The most important aim of heavy ion experiments performed at the RHIC accelerator 
besides the production of quark matter is the reconstruction of its properties from the 
observed hadron measurables. 

The investigations of this problem are focused around three main questions: a) 
the effective mass of the quarks, b) the momentum distribution of quarks, c) the 
hadronization mechanism. 

Regarding the quark mass in quark matter up to now essentially two different 
assumptions has been made: i) the quark matter consists of an ideal gas of non- 
interacting massless quarks and gluons and this type of quark matter undergoes a first 
order phase transition during the hadronization; or ii) the effective degrees of freedom 
in quark matter are dressed quarks with an effective mass of about 0.3 GeV, which 
mass is the result of the interactions within the quark matter; this type of quark matter 
clusterizes and forms hadrons smoothly, following a cross-over type deconfinement - 
confinement transition. 

Considering the momentum distribution in quark matter one meets the following 
ideas: The bulk part of the quarks is produced in soft processes, these quarks have a 
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fairly exponential transverse momentum distribution on the top of a collective transverse 
flow in heavy ion collisions. On the other hand, quarks with large transverse momentum 
(pr > 2 — 4 GeV) are viewed as produced mainly in hard collisions and their spectra 
obey a power law shape described by the functional form of fragmentation functions. 

Finally one meets the worthy problem of hadronization. In early publications QjJ E] 
the hadronization was assumed to happen through a quasi-static, first order phase 
transition. However, as the experiments have suggested that the whole lifetime of the 
heavy ion collision is very short, other non-equilibrium dynamical models obtained more 
and more credit. One of the most successful description is based on quark coalescence 
[31 IU El EH HI IE] • This idea is supported by recent RHIC experiments 

In this paper we present a new model of the late phase of the heavy ion collision. In 
this approach the effective quarks are massive and have a particular mass distribution 
described in Section 2. In Section 3 the quark momentum distribution is studied in 
details. Section 4 summarizes the appropriate expressions for the coalescence processes 
in this environment. Numerical results and comparison to experimental data are 
displayed in Section 5. We discuss the obtained results in Section 6. 

2. Mass distribution of Quarks 

We assume that in hot and dense deconfined matter the average interaction of the 
massless (current) quarks generate an effective mass for them. Several analyses of lattice 
QCD calculations supports the appearance of massive quark degrees of freedom [TUIITT] . 
(According to these analyses the gluons become very heavy just before hadronization, 
and thus we neglect their explicit degrees of freedom in the models.) The mass of these 
quasi-particles is in the range of 0.3 — 0.4 GeV around the phase transition temperature. 
In the early approaches the quasi-quarks have a sharp, momentum independent mass. In 
a realistic approach, however, there is no reason to assume that the interacting quarks 
behave as on-shell particles. Due to the residual interaction the mass of the quasi- 
quarks may be distributed in an interval of about a few hundred MeV while having a 
maximum in the range of 0.3 — 0.4 GeV. The shape of this distribution is presently not 
known f. One may want to investigate the effect of different spectral functions on the 
hadronization process. In the present work we consider the following mass distribution: 



This distribution vanishes both for too small and too large masses and has a smooth 
maximum around the familiar valence quark mass value (//). FigQ displays this mass 
distribution for /i = 0.25 GeV and T c = 0.26 GeV. 

| Not thinking that quarks are dressed only by one partner particle (gluon), there is also no reason to 
stick with the Breit-Wigner form for this distribution. 
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Figure 1. Distribution of the effective quark mass in eq. 



3. Momentum distribution of quarks 

In the RHIC experiments it was observed, that the produced particles show a nearly 
exponential spectrum at low transverse momenta and a power-law like one at high 
transverse momenta. This distribution was considered in earlier publications as a sum 
of three components produced by different physical mechanisms: a) thermal equilibrium, 
b) coalescence, c) fragmentation. The relative weights of these components were 
determined by a fit to the experimental data. In the present work we assume an analytic, 
parameterized transverse momentum distribution, which is similar to an exponential 
distribution for low momenta and a power law distribution for large momenta. The 
phase space density of effective quarks with a variable mass m, transverse momentum 
PTi coming from an angular direction <fr and a coordinate rapidity rj is given by 



Here q is a parameter characterizing a possible deviation from the usual Boltzmann 
distribution, which is recovered for q = 1. The parameter T is in general not the usual 
temperature, but in the Boltzmann limit, q — > 1 it becomes the canonical temperature. 
That which values of q and T fit experimental heavy ion data best, may depend on the 
transverse momentum range considered: a q value not equal one fits usually to much 
higher momenta than the original Boltzmannian assumption. 

Although the expression in eq. (J2J) coincides with the form of the Tsallis 
distribution jT^j (as many simple minded fits do), we do not necessarily have to assume 
a special equilibrium state, which is considered in the non-extensive thermodynamics. 
One may assume such an equilibrium Tsallis system in quark matter if one of the 
following conditions is fulfilled: i) there is an energy dependent random noise in the 
system not reflecting a constant temperature; ii) the system contains a finite number of 
particles which interact indefinitely long; iii) there are long range forces in the system. 
We realize, that such conditions may be found in a fireball produced in high energy 
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heavy ion collisions, it is just not sure, whether the time is enough to arrive at such a 
special equilibrium state. 

The energy E, appearing in eq.(J2J) is given by the blast wave formula: 

E(p±, m, 0, 77) = i(\jm 2 +p\ ■ cosh(r/) - v fiow ■ p± ■ cos(»). (3) 

assuming a homogeneous transverse flow velocity, Vfi ow , with a corresponding Lorentz 
factor, 7. In eq.Q the effect of an average "blue shift", caused by the hydrodynamical 
flow, is taken into account. 



4. Coalescence of Massive Quarks 



Using the covariant coalescence model of Dover et al. |T3j , the spectra of hadrons formed 
from the coalescence of quark and antiquarks can be written as 

dN h dN h 



E ■ 



9h 



(Ph ■ da h:ll )F h (x h ;p h ). 



(4) 



d 3 p dyp T dpTd(J)p (2tt) 3 

Here Fh(xh',Ph) is an in principle eight dimensional distribution (Wigner function) 
of the formed hadron. Furthermore da denotes an infinitesimal element of the space- 
like hypersurface where the hadrons are formed, while gh is the combinatorial factor for 
forming a colorless hadron from a spin 1/2 color triplet quark and antiquark. For pions 
the statistical is g n = 1/36 , for antiprotons g-p = 1/108 [6 . 

Assuming an instant hadronization at a sharp proper time, r, with cylindrical 
symmetry one gets 

- — = j^-t f t ■ m±cosh(y - r])dr]rdrd(f)F h (x h ;p h ) 
dypxdpTdcpp (zTry J 

Here y is the momentum rapidity of the hadronic four-momentum p^ and 77 is the 
coordinate-rapidity of the formation point. For y = and = eq.(JH) becomes 



dN h 



dyp T dp T c 



9h 



rA v /2 

tiR 2 t -m ± - dr] cosh(?7) F h (p h ; 77) (5) 

JA v /2 



y=o ( 2 -) 

Assuming that we get contribution to the spectra mainly from the << 1 region 
eq.(jnj) simplifies to 

dN h 



dyp T dp 7 



y=0 



_9h_ 

(27T) 



ttR 2 tA v ■ m± ■ F h ( Ph ; 0) 



(6) 



where m±_ = y m 2 + p\ is the transverse mass, (f) p the angle of the transverse momentum 
of the hadron and R is the cylindrical radius. Introducing the hadronization volume, 
V = A v tR 2 ti, and the abbreviation Fh(pj_, 0) = Fh(p±) we finally get 

dN h 



dyp T dp q 



y=0 



9h 

(27T) 



■Vm ± F h (pj 



(7) 
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Assuming that the mesons are produced by coalescence of particle a with mass m a 
and particle b with mass m& , the meson source function, F/v/(p M ), will have the form: 

F M (Ph) = J d 3 p a d 3 p b f a (p a ;0)f b (p b ;0)C M (Pa,PbiPM)- (8) 

In eq. (JHJ) the coalescence function CuijPaiPbiPM) is the probability for particles with 
four momenta p a ,Pb to form a meson with four momentum pu- We assume a Gaussian 
wave-packet shape being maximal at zero relative momentum (characteristic for s-wave 
internal hadronic wave functions): 

C M (Pa,Pb,PM) = % ■ e-^- p ^ 2 (9) 

The parameters ctu arid P c reflect properties of the hadronic wave function in 
the momentum representation convoluted with the formation matrix element. For the 
present consideration we assume that the system consist of a mixture of effective particles 
having different masses. The densities of particles with mass m are weighted with the 
corresponding mass distribution function, p(m), eq.(JT]). 

First we deal with the formation of mesons. We assume that the total mass and 
momentum of a hadron will be combined additively from that of the constituents, but 
P c is so small, that that practically zero relative momentum partons form a meson. This 
leaves us with 

Pa = Pb = Pm/2 (10) 
m a + m b = m M (11) 

This leads to the coalescence function 

C M = olm -5(p a - Pm/2) -S(p b - p M /2) ■ 5(m a + m b - m M ) (12) 

. Thus we arrive at the following meson distribution function: 

rm M rm M 

Fm{Pu 0) = a M • / dm a I dm b ■ 8{m M ~ (jn a + m b )) 
Jo Jo 

■ pijn a ) ■ f q (p t /2;m a ) 

■ p(m b ) ■ f b (p t /2;m b ). (13) 
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Figure 2. The product of the two quark-mass distributions, p{m) ■ pira-^ — to), in 
ea l|13|) . The maximum of the product is at to = 
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Fig|21 shows the convolution of two quark-mass distribution functions appearing 
in the coalescence probability eq (jl3|) for the coalescence production of a pion. 

One observes, that the product of the two mass distributions always has a maximum 
around m = Already the naive coalescence approach works, whenever this 

maximum is sharp enough. 

Similar argumentation leads to the three-fold coalescence expression. In the place 
of eq.(JTTJ) we get: 

Pi = P2 = Ps = Pb/3 (14) 
mi + 1^2 + ™3 = m B (15) 
and the baryon source function becomes 

rmp r rm pr rm pr 

F~B{Pt) = ol b - dm 1 / dm 2 / dm 3 
Jo Jo Jo 

■ pN/ 9 (PtA m i)) 
• P(m 2 ) f q (p t /3,m 2 )) 

■ p(m 3 ) f q (p t /3,m 3 ))5(m B - (m 1 + m 2 + m 3 )). (16) 




Figure 3. Product of three quark-mass distributions, p{mi) ■ p(m 2 ) ■ p(m 3 ), from 
eq. <|16l) . Horizontal axes connect to masses mi and ni2. Mass 7713 is defined as 
m 3 = nip — {nil +ni 2 ). 



Fig 01 shows the product of the three quark- mass distribution functions appearing 
in eq.(HI)J). 
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Figures 2 and 3 shows that the maximum contribution to the coalescence integrals 
are obtained from the equal mass part of the mass distributions. 

At the end of this section we may compare the basic philosophy of the different 
coalescence models. In all cases we met the problem of mismatch of initial and final 
state quantum numbers. To resolve this problem different solutions are proposed: 

a) In the quantum mechanical approach the overlap of initial and final state 
wave function is calculated. In this process - simply said - the quantum mechanical 
uncertainty bridges the mismatch gap. 

b) In the most often used classical approximation the coalescence probability is 
smeared around the momentum conservation, Ref. [T3*] . 

c) In the present approximation in the coalescence transition the momentum is 
strictly conserved, as in Ref. [Zj, but the masses of the particles of the initial state have 
a finite width distribution. 

5. Numerical results and comparison with experimental data 

The parameters used in the calculations are as follows. \i = 0.25 GeV, T c = 0.26 GeV, 
T = 0.07 GeV, q = 1.19, Vfi ow = 0.56. The % and as coalescence constants were 
adjusted to fit the absolute value of the measured 7r and p yields. In Fig.(HJ) both the 
7T~ and 7T° data are displayed. 

In the Figures 4-6 we show the calculated distributions together with the 
corresponding experimental data. The experimental points are taken from Ref. [T%] . 
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Figure 4. Measured pion spectra in central Au+Au collisions at y/s = 200 AGeV Jl] . 
Full line indicates the calculated pion yield. 
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Figure 5. Transverse momentum spectrum for antiproton production in central 
Au+Au collisions at y/s — 200 AGeV Full line indicates our calculation. 
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Figure 6. The p/ii~ ratio as a function of pr on a linear scale in central Au+Au 
collisions at y/s = 200 AGeV [Hj. 
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6. Conclusion 

In earlier publications [BJII] it was assumed that for the low, medium and high transverse 
momentum intervals the hadrons are produced with different mechanisms, via thermal, 
coalescence and fragmentation processes. In the present work we proposed a model, in 
which, besides the momentum, the energy is also conserved to a good approximation 
in the coalescence process. This model yields results, which agree fairly well with 
experimental heavy ion data in the transverse momentum range from about 1 GeV to 8 
GeV, (see Figs. (4) - (6)). This agreement supports the basic assumptions of our model: 
a) the masses of the partons have a finite width distribution (Fig.l); b) the momentum 
distribution of partons is described by one and the same function for a wide range of 
transverse momenta; c) the hadronization happens via effective quark coalescence; d) 
in the hadronization process both momentum and energy are nearly conserved locally. 

Finally we note that a spherical inhomogeneity in the momentum distribution of 
the quarks leads to an azimuthal correlation of hadrons, if they are produced by the 
coalescence mechanism described in the present paper. This effect may be worth of a 
closer look in future studies. 
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